Quantum optics in different representations of the algebra of canonical commutation 
relations (I): Unexpected properties of Rabi oscillations — theory and experiment 
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We discuss the Jaynes-Cummings model in different representations of the algebra of canonical 
commutation relations. The first conclusion is that all the irreducible representations lead to equiv- 
alent physical predictions. However, the reducible representation recently introduced as a candidate 
for 'QED without infinities' leads to new effects. We analyze from this perspective the experiments 
on Rabi oscillations performed by the Kastler-Brossel Laboratory group from Paris. Surprisingly, 
the results seem to support the reducible representation approach. We also discuss possibilities of 
more definitive tests of the new formalism. 



I. INTRODUCTION 

It is known that different representations of the same 
Lie algebra correspond to different types of physical 
systems. Parameters that characterize the representa- 
tions have a meaning of quantum numbers (angular mo- 
mentum for the rotation group, mass and spin for the 
Poincare group, and so on). In principle, if we find a new 
representation of some physically meaningful group or al- 
gebra we should seriously consider the possibility that the 
representation corresponds to a physical system, perhaps 
yet unknown. This is how the anyons were predicted [l| . 

In quantum optics the central role is played by the Lie 
algebra of canonical commutation relations (CCR) 



[a(k),a(k')*\ = 5{k,k')I(k) 



(1) 



The element I(k) commutes with all the elements of 
the algebra and, by the Schur lemma, is proportional to 
the identity operator if the representation is irreducible. 
Fields with different boundary conditions correspond to 
different I(k) and S(k,k'), and thus to different repre- 
sentations of CCR. 

Still, the problem seems deeper than that. First of 
all, the theorem of von Neumann states that there 
exists an infinite number of incquivalcnt irreducible rep- 
resentations of CCR for systems with infinite numbers 
of degrees of freedom (here corresponding to the infinite 
number of different wave vectors). It is not clear which 
representation to choose, even if one keeps boundary con- 
ditions fixed and restricts the analysis to free fields. Sec- 
ondly, as stressed by Dirac in his last published paper 
0, one should not ignore the physical potential inher- 
ent in reducible representations. These are unavoidable 
in quantum physics since tensor products of irreducible 
representations are not themselves irreducible. Compos- 
ite quantum systems are inherently reducible. 

One of the representations that gained particular pop- 
ularity in quantum optics employs the Hilbert space of 
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infinitely many harmonic oscillators. To each frequency 
of the field there corresponds a separate oscillator. Phys- 
ically the representation is rather pathological, just to 
mention the problem with the infinite energy of the 
ground state. Mathematically, the representation is also 
pathological in the sense of involving a non-separable 
Hilbert space typical of infinite tensor products Q, a 
fact rarely mentioned in the quantum optics literature. 

It is relevant to mention in this context that recently, 
in a series of papers 0, IE S 13 one of us investigated the 
possibility of electromagnetic field quantization in terms 
of certain reducible representations of CCR. The main 
feature of the proposed formalism was that there was 
no link between the number of frequencies allowed by 
field boundary conditions and the number N of oscilla- 
tors used in construction of the Hilbert space. In par- 
ticular, it made perfect sense to speak of fields modelled 
by a finite number of oscillators even if the number of 
frequencies in the field was infinite. And vice versa, a 
monochromatic field could be modelled with any N. 

It is easy to understand why this is possible if one 
thinks of each oscillator as a wave packet involving all 
the possible frequencies. The parameter iV occurring in 
the representation has a status of a quantum number, 
and the vacuum is a Bose-Einstein condensate of N os- 
cillators at zero temperature. The limit N — > oo plays 
a role of a correspondence principle mapping the new 
formalism into the standard one. What is interesting, 
when it comes to computing averages it turns out that 
the wave function of the ground state plays a role of a 
cut-off function regularizing integrals at ultraviolet and 
infrared regimes. The new formalism thus automatically 
introduces many elements that are imposed in an ad hoc 
manner in the standard one. 

However, if in reality N is large but finite then quan- 
tum field theory we know nowadays is the N = oo ap- 
proximation of a more fundamental, and hopefully more 
consistent, N < oo theory. The question is: Why the 
N = oo approximation works so well, and where to look 
for experimental manifestations of a finite iV? 

In this paper we give partial answers to both parts of 
the question. We concentrate on the Jaynes-Cummings 
model where we solve Heisenberg equations of motion for 
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the atomic inversion operator. A difference with respect 
to the usual treatment of the problem is that we begin 
with a representation-independent level. The first con- 
clusion we find at that stage is that all the irreducible 
representations of produce physically equivalent re- 
sults. 

As the next step we choose the reducible representa- 
tion parametrized by N < oo and compute evolution 
of atomic inversion with different vacuum, thermal, and 
coherent-state initial conditions. 

Here the situation changes. But first of all it has to be 
stressed that in the limit TV — > oo we recover the stan- 
dard formulas. This fact is of crucial importance for the 
whole approach since it plays a role of a correspondence 
principle and guarantees that the new theory may be re- 
garded as a generalization of the one based on irreducible 
representations. 

For finite N there are differences which we compare 
with experimental data of 9]. What is interesting, it 
seems that the data are more consistent with a finite iV 
than N = oo (N = 280 for the nearly vacuum state at 
T = 0.8 K, N = 420 for the coherent state with 0.4 
and 0.85 photons in average). The main reason why the 
experiment does not fully support the standard theory is 
that for finite N one expects a faster relaxation of Rabi 
oscillation than what one might expect on the basis of 
N — oo, and this is precisely what seems to happen in 
the experiment. 

The relaxation occurring for finite N is not a decay but 
a beat and therefore waiting sufficiently long one should 
see a revival of the Rabi oscillation. We show that cavi- 
ties with lifetimes of a few hundred /is in principle allow 
for observation of the revival. Of particular interest are 
the maser and mazer experimental setups [Tol UlTIT^. Il3| . 
but we leave it for a future work. 

If our intuitions are correct and the finite N represen- 
tations are more physical than the limiting TV = oo case, 
then many additional questions have yet to be answered. 
In the paper that accompanies the present one we address 
the issues of Lorentz and gauge covariance and test the 
formalism on another exactly solvable model: Quantum 
fields produced by a classical current. 

The present paper is organized as follows. In Sec. II 
we discuss the Jaynes-Cummings model at a representa- 
tion independent level. In Sec. Ill we assume that the 
representation is irreducible and conclude that all such 
representations imply equivalent physics. In Sec. IV we 
switch to the reducible representation involving N < oo 
oscillators. Then, in Sec. V, we write the explicit form of 
the inversion operator Rs(t) in this representation and, 
in Sec. VI, compute atomic inversion with different ini- 
tial conditions. In particular, we discuss in detail the 
limit N — > oo and its links to the law of large numbers 
and renormalized parameters. Some technicalities and 
remarks about the reducible representation are moved 
into Appendices. 



II. REPRESENTATION INDEPENDENT 
FORMULATION OF THE JAYNES-CUMMINGS 
MODEL 

The Jaynes-Cummings model represents a two- 

level atom interacting with a single mode of electromag- 
netic field in a cavity. The cavity boundary conditions 
imply that the set of free-field momenta is discrete and 
it is convenient to work from the outset with the discrete 
notation 



[aj,a%] = 5 jk Ij. 



(2) 



Here 8jk is the Kronecker delta, the asterisk denotes the 
adjoint in the sense of *-algebras, and I* = Ij. We as- 
sume there exists a free-field Hamiltonian Hq satisfying 



[aj,H ] 
[a*, H ] 



(3) 
(4) 



Note that Hq cannot, in general, be given by 
the latter works only if Ij is an identity, which we do not 
assume at the present stage. 

Let us now select a frequency Uj = lu and assume 
that only this frequency couples to the two-level system. 
The corresponding CCR operators will be indexed by lu, 



that is: 



Jo 



a u , a 



Jo 



* T 



We also split 



Hq into two parts: Hq commuting with and a* , and 



Hq = luN^, where 



(5) 
(6) 



Alternatively, we may assume that iV w exists and define 

Hq - LuN u . 

The model is given by the full Hamiltonian 



H = w R 3 + Hq + gR+a u + gR-a* 



(7) 



We employ the usual notation [15J where Ri = ai/2, 
R± = Ri ± ii?2, c; are the Pauli matrices, and g is a 
complex coupling parameter. It is useful to split H into 
three mutually commuting parts: 



H = Hq + lojV + n, 

O = AR 3 + gR + a u + gR-a* 



(8) 
(9) 
(10) 



A = luq — lu is the detuning. The evolution operator 
factorizes: 



Vi 



-iClt 



cos(il/{f) — i 



sin(f^i) 
VP 



n. 



X = RHuj + a* a w + 



(11) 
(12) 

(13) 
(14) 



3 



Let us note that the operators M — R 3 I U + a* a w , oc- 
curring in X, and j\f = R$ + N u , occurring in H, should 
not be identified, although both are equivalent in irre- 
ducible representations of CCR, and both commute with 
H independently of the choice of representation. Con- 
structing M. and M in the reducible representation we 
shall explicitly show that they are essentially different. 



Only (|12() is relevant in the context of Heisenberg- 
picture dynamics of R 3 (t) since [R 3 ,Hq] = [R 3 ,J\f] = 0. 
One can verify by a straightforward calculation that 
idUt/dt — UtH, and it is clear that the solution is valid 
for any representation of the Lie algebra (0) . Employing 
(11211 one finds the Heisenberg picture evolution of R 3 : 



R 3 (t) = R 3 [l-2\g\ 2 X 



.sin 2 (tt R t) 



n 2 



A sin 2 (Q R t) 



.sin(2fi fl i)\ p 



Asin 2 (0^) .sm(2tt R t) 



2n, 



(15) 



The proof of (|15JI is outlined in the Appendix. 

The next important notion that can be introduced at 
a representation independent level is the displacement 
operator. The operator is defined in the usual way as 

D{z) = exp^2(z j a* j - Zjttj) (16) 
3 

and satisfies 

D{z)*a 3 D{z) = aj +ZiIi, (17) 

D(z)*a*D(z) = aj + Zjlj, (18) 

D(zyijD{z) = Ij. (19) 

Acting with D(z) on a vacuum vector we obtain a coher- 
ent state. Its form depends on what is meant by vacuum 
in a given representation. 

III. EVOLUTION OF ATOMIC INVERSION 
OPERATOR IN IRREDUCIBLE 
REPRESENTATIONS 

Assume we work in an irreducible representation with 
some carrier Hilbert space containing a vector |0) annihi- 



lated by all a u . The abstract ^-conjugation in the algebra 
can be replaced by Hermitian conjugation of operators in 
the representation. The representation satisfies 



[a j ,al]=ZS jk l, (20) 



where Z is a real positive number (by Schur's lemma an 
element that commutes with all elements of the algebra 
is a constant times identity if the representation is irre- 
ducible, i.e. Ij = Zl for some Z, and for all j; had Z 
been negative we would have called atj a creation opera- 
tor, taken vacuum annihilated by at, and appropriately 
adjusted the notation). The solution (|15fl involves the 
operator \g\ 2 X whose representation reads 

| 5 | 2 A = \g\ 2 (ZR 3 + ata UJ + Z/2) (21) 
= |.g| 2 (i?3 + 4^ + 1/2) = \g\ 2 X, (22) 

where a u — a^/y/Z, [a,j, &£] = Sjkl, and g — \fZg. Now, 



R 3 (t) = R 3 \\-2\g\ l X 



- sin 2 (ft fl .f) 



n 2 



A sm 2 (n R t) 



fl 2 



.sin(2fi fl f)\ . 



/\sin 2 (Q R t) sixi(2Q R t)\- 

+ i — — fgR-a^, 



n 2 



2flj 



(23) 



where f2 = y A 2 /4 + |g| 2 X. Clearly, the only difference 
between different representations is in the values of the 
coupling constant g = g\f~Z. Since g is proportional to 
the electron charge, different representations effectively 
differ by the value of the electron charge eo = eoV~Z. 



This type of rescaling is exactly what occurs in transi- 
tion from the bare charge eo to the physical, renormalized 
charge e p h = eo \fZ3~ typical of renormalized electromag- 
netic fields. Let us note finally that 

5jZ)(z)|0) = ZjVZD(z)\0) (24) 
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and thus the relation between the parameter z and the 
'physical' amplitude z is also renormalized: z — z\[Z. 
After these rescalings of the 'bare' parameters we will 
obtain identical formulas for coherent state averages of 
i?3(t), independently of our choice of the irreducible rep- 
resentation. We conclude that the theorem of von Neu- 
mann does not bring to our problem anything physically 
important. 

So let us switch to reducible representations. 



IV. TV < co REDUCIBLE REPRESENTATION 

The representation is constructed as follows. For sim- 
plicity we ignore the polarization degree of freedom. Take 
an operator a satisfying [a, a^] — 1 and the kets \k) cor- 
responding to standing waves in some cavity. We define 



<z(fc) = |fc)(fc| Oa, I(k) = \k)(k\ <g> 1. 



(25) 



The operators Il25[i satisfy Q), where S(k,k') is the 3D 
Kronecker delta. The fact that I(k) is not proportional 
to the identity means that the representation is reducible. 
In our terminology this is the 'TV = 1 representation'. Its 
Hilbert space Ti is spanned by the kets \k,n) = \k)\n), 
where a)a\n) = n\n). Such a Hilbert space represents 
essentially a single harmonic oscillator of indefinite fre- 
quency [5| . An important property of the representation 
is that J2k ^CO = ^ is the identity operator in Ti. A 
vacuum of this representation is given by any state anni- 
hilated by all a(k). The vacuum state is not unique and 
belongs to the subspace spanned by |fe,0). In our nota- 
tion a N = 1 vacuum state reads \0) = J2 k O(k)\k,0) 
and is normalized by |0(fc)| 2 = ^ fc Z(k) = 1. Such a 
vacuum is exactly analogous to a single-oscillatorground- 
state center-of-mass wavepacket. As shown in ja,l3 (£° r 
a more complete discussion see the companion paper (8j) 
in a fully relativistic formulation the maximal probabil- 
ity Z — maxfc{Z(fc)} is a Poincare invariant, an aspect 
worth keeping in mind. For TV > 1 the representation 
space is given by the tensor power Tt = Ti® N , i.e. we 
take the Hilbert space of TV (bosonic) harmonic oscilla- 
tors. Let A : TL — > Ti be any operator for TV = 1. We 
denote = I®^- 1 ) ® A® J®^-*), A (n ) : H -> H, 
for 1 < n < TV. For arbitrary TV the representation is 
defined by 

1 N 

^ = ^E a ( fc ) (n) ' ( 26 ) 



n—1 

i N 

n=l 

[a(k),a{k')^} = 5(k,k')l(k), 

E^( fc ) = L = i 0N 



(27) 

(28) 
(29) 



and the TV-oscillator vacuum is the TV-fold tensor power 
of the N = 1 case, a kind of Bose-Einstein condensate 



consisting of TV wavepackets: 

\0) = \0)®---®\0) = \0)® N . 



(30) 



The free- field Hamiltonian is, for N—1 and u>k = \k\ 



H o = Efe w fe a(fc) t a(fc) = 



,t, 



(31) 



In each eigensubspace with fixed |fc) the operator H is 
just an ordinary Hamiltonian of the oscillator with fre- 
quency uJk- For arbitrary TV the generator of free field 
evolution is the Hamiltonian of TV noninteracting oscil- 



lators, i.e. H_ Q = J2n=i^a l> - Let us stress that Hq 
should not be confused with J2k Wfea(fc)' l 'a(fe) which does 
not have a clear interpretation in this context and does 
not describe noninteracting oscillators. We have seen al- 
ready that the operator a(kya(k) will nevertheless play 
an important role in the Jaynes-Cummings problem. Our 
definition of iJ implies that [a(fc), H Q ] — tOka(k) which 
is the formula we required at the representation inde- 
pendent level. Let us now select some p with frequency 
u> = \p\ to be the mode that interacts with the two- 
level atom. To simplify further calculations we denote 

a{ P ) = a,,, Up) = 4,, \p)(p\ ®i = Pi,p = i-Pv 

The projector Pj is related to by 



(n) 



N\ 



I- 



I 



1 1 



One recognizes in (|32|l the frequency-of-success operator 
discussed in the context of quantum laws of large num- 
bers [3 EHIl E IM IHHll. Eigenvalues of 1^ coincide 
with all the possible frequences of 'heads' in TV trials of 
coin tossing, i.e. 1^ = J2^=o( s /N)P(s/N) with spectral 
projectors 



P(s/TV) 



E 

S1 + ... + SJV— s ': s j— 0,1 



(33) 



The explicit form l|33|) shows that P(s/N) commute with 
all a(k), a(ky , and I_(k). Of particular importance is the 
splitting of the Jaynes-Cummings problem into subspaces 
corresponding to a given s/N. We define 

a(s) = a^P(s/N), (34) 

a(s) f = 0lP(s/N), (35) 

and obtain the representation 

[a(s),a{s^} = (s/N)P(s/N) (36) 

of CCR in the Hilbert space H(s) = P{s/N)H. 

Each TL(s) is an invariant subspace for the Jaynes- 
Cummings dynamics, and in each such a subspace we 
effectively deal with a representation given by i|2U[l whose 
Z = s/N, as we shall see in the next Section. 



V. EVOLUTION OF ATOMIC INVERSION 
OPERATOR IN TV < oo REPRESENTATION 

As the general formula l|15fl is valid also in this rep- 
resentation and all the operators that occur there com- 
mute with P(s/N), we begin with splitting Pa(t) into 
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parts acting in the invariant subspaces TC(s). We employ 
the usual notation R 3 = (|+)(+| — |— )(— 1)/2 and treat 
|+)(+| + | — )(— | as the identity operator (to be more ex- 
act we should tensor the atomic-space identity with the 
field-space identity, but we prefer this simplified conven- 
tion) . 

Denoting X(s) = XP{s/N) we get X = J2s=o x ( s ) 
and 

X(s) = (±R 3 + a( S ?a( S ) + ^)p(s/N) (37) 
= ^\+){+\P(s/N) + a(s^a(s). (38) 

Commutation relation (|36(1 implies that eigenvalues of 
a(sya(s) are sn/N, n = 0, 1, 2, . . . . In spectral represen- 
tation 

oo 
n=0 

The spectral projectors satisfy H(n, s)P(s/N) — n(n, s), 



Er=o^( n,s ) = P( s /N). Spectral representation of 
X(s) therefore reads 

OC 

X ^ = ^E(i+)<+i+ n ) n (^ s ) ( 4 °) 

oo 

= ^E nfl K s ). (41) 

Spectral projectors tl(n, s) of X(s) are related to II(n, s) 
of a(sya(s) by 

fl(n,s) = |+)(+|n(n-l, S ) + |-)(-|n(n,s), (42) 

for n > 0, and n(0, s) = |-)(-|II(0, s). In the Appendix 
we show that [R + a(s),fl(n, s)] = [i?_a(s)t,n(n, s)]. 
Evolution of the atomic inversion operator is given in 
our reducible representation by 



N oo 



R 3 (t) = R 3 -Y,Y, 2R ^\ 



s=0 n=0 
N oo 



, ns sin 2 (ty/A 2 /4 + \g\ 2 ns/N) 
iV A 2 /4+ \g\ 2 ns/N 



U(n,s) 



^ ^ , A sin 2 (VA 2 /4 + \9?n S /N) .sin (2*^/4 + \g\ 2 ns/N) \ 
l^Lu I I u 2 _.;»r 1 » . 1 w rrr Mn,s)gR+a{s) 



s=0 n=0 
N oo 

EE 



A 2 /4 + \g\ 2 ns/N 



2^A 2 /4+ |5| 2 ns/iV 



A sin 2 (VA 2 /4+ | 5 | 2 ns/iV) . sin (2^A 2 /4 + \g\ 2 ns/N) 



A 2 /4+\g\ 2 ns/N 



+ l- 



2 v /A 2 /4+ \g\ 2 ns/N 



gR-a{s) f fl(n,s). (43) 



Formula (|43J) can be directly compared to 12.'{[l by means 
of the following simple rule: Skip the sum over s and 
set s/N — 1. An alternative recipe is to skip the sum 
over s and set s/N — Z, which corresponds to CCR 
with 21 at the right-hand-side. Yet another intuitive 
rule can be found by means of the law of large numbers 
and works for the weak limit N — > oo. To make it precise 
one has to switch to the level of averages. We will see 
that the weak law of large numbers plays a role of a 
correspondence principle between our N < oo formalism 
and the standard regularized one. 

VI. EVOLUTION OF ATOMIC INVERSION IN 
TV < oo REPRESENTATION 

Acting on the vacuum state l|3()|l with the displacement 
operator 1)16(1 one obtains a coherent state. Here we are 
interested in a monochromatic coherent state with fre- 
quency u> 

k) = exp(zal-zg^)\0). (44) 



The coherent state is not an eigenstate of the annihilation 
operator but a direct sum of its eigenstates. Indeed, 
the decomposition 

N N 

\z) = ^2P(s/N)\z)=^2\z(s)) (45) 

s=0 s=0 

accompanying ^ = J2 s =o^ P ( s / N ) = Eflo a ( s ) in> 
plies that 

a(s)\z(s)) = (s/N)z\z(s)). (46) 

The analogy of the latter eigenvalue problem to (|24fl is 
evident. Alternatively, one can say that the coherent 
state is a generalized eigenvector of a^, i.e. 

sj*) = *L|2>- ( 47 ) 

Another state of interest, particularly in the context of 
experiments, is the mixture 

P = (p+l+X+l+P-l-X-O/Ofieid- (48) 
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A. Reduced inversion operator: Coherent state 



involving only the atomic degrees of freedom: 



Evaluating an average of (|43|l in an arbitrary coherent 
state \z) wc obtain the reduced operator 



R z (t) = (z\R 3 (t)\z) 



(49) 



N 



R z {t) = i? 3 + ^-^(i|n(v)k: 



N 

71=0 



x -\g\\n + l)- 



{t^A 2 /A+\g\ 2 {n + l)s/N) 
A 2 /4+ \g\ 2 {n+ l)s/N 



+X+l + l5| 2 » 



sin 2 (i^/A 2 /4 + \g\ 2 ns/N) 
A 2 /4+ \g\ 2 ns/N 



l-X-l 



zg 



( A sin 2 (t ^A 2 /4 + |.g| 2 (n + l)s/N) . sin (2t ^/MJi+ \g\ 2 (n + Tjs/N) 



A 2 /4+ \g\ 2 {n + l)s/N 



2^A 2 /4 + \g\ 2 (n + Tjs/N 



R4 



_ ' A sin 2 (t v / A 2 /4+ \g\ 2 (n + Tjs/N) sin (2i v /A 2 /4 + \g\ 2 {n + Tjs/N) ' 
+ zg I — ao i a . i —, — . — —r 1- i = — ]R- 



A 2 /4+ \g\ 2 {n + l)s/N 



2y/A 2 /4:+\g\ 2 (n + l)s/N 



(50) 



where (see Appendix) 
(z\U(n,s)\z) 



(51) 



-\zJs/N\ 



N 



N- 



Here 



(o\UQ) = (o\i u \o) = \o{ P f 



(52) 



is the probability of finding the momentum p correspond- 
ing to the resonant mode u> = \p\, if the vacuum state of 
the field is \Q). In the last term of 151fl one recognizes the 
binomial distribution for N Bernoulli trials with proba- 
bility of success equal to Z^. 



B. Vacuum-state initial condition: jz) = \Q) 

Assume initially there are no photons and the atom is 
in either ground or excited state, i.e. 



m = \±)\0. 



(53) 



In this case only the n = term counts in l|5U|) . The case 
of atomic ground state, \^f) = |— )\0), is trivial since 

w(t) = (*|ife(*)l*> = HRz(t)\-) = -1/2. (54) 

However, starting with the excited state \*f>) = |+)|0) wc 
find 



JV 



w(t) = (ms(t)m = (+\R*(t)\+) = o - E \af 



v / A 2 /4+\g\ 2 s/Nt ( N 



s=0 



N 



A 2 /4+\g\ 2 s/N 



- ZJ) 



N-s 



(55) 



r 



The law of large numbers for the binomial distribution 
implies that 



r M 1 | ,2 7 sin 2 y /A 2 /4+\g\ 2 Z LJ t 



A 2 /4+| 5 | 2 Z a 



i.e. the frequency s/N approaches the probability of suc- 
cess in a single trial, that is s/N — > Z u . 

It is instructive to compare this result with the one 
we would have found had we started with the general 



irreducible representation whose right-hand-side is Zl. 
The corresponding result reads 



w(t) 



„ sin 2 ^JA 2 /±+\g\ 2 Zt 



\9\*Z 



A 2 /4+\g\ 2 Z 



(57) 



These two formulas appear so similar that one can easily 
overlook an important difference: The parameter Z in 
(|57|l is a number independent of lo, whereas Z u in i|56|) 
is a value of a function vanishing for lu — > oo. Charge 
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renormalization cannot be performed in the same way in 
both models, because charge must remain a relativistic 
invariant. 

In the irreducible case the recipe is simple: e p h = 
egy/Z. In the reducible case one first extracts from 
the relativistic invariant Z — maxfe{Z(fc)} 0,0, 0- Then 
one writes Z u = Zxw, redefines charge e p h = eov^Z, and 
finally 



lim 



»(*) = \ - 



l5 P h| 2 Xc 



sinVA 2 /4+|g P h| 2 X^ 
' AV4 + |g ph | 2 ^ 



(58) 



The solutions l|57|l and (|58|l are identical up to the pres- 
ence of the cut-off function x^ m • Needless to say 
the function would necessarily regularize the interaction 
if we decided to work with extremely high frequencies u>. 



Let us note that we did not introduce any cut-off in 
the Hamiltonian. The cut-off has appeared automati- 
cally through the structure of the vacuum state typical 
of the reducible representation. The theory gets regu- 
larized even though we do not really need it in such a 
simple example. The regularization neither solves here 
any problem nor spoils anything. If we assume that for 
optical frequencies Xui ~ 1 we obtain exact agreement 
between the irreducible case and the N — > oo limit of the 
reducible one. This is an example of the correspondence 
principle we have mentioned in the introduction. 

The next step is to understand if one really needs N = 
oo, and if some finite N cannot, in fact, be consistent 
with experimental data. We first rewrite l|55|l by means 
of the renormalized coupling 



w(t) 



" | 3ph | 2 s sin 2 v AV4 //pi: VlZ.Y)/ ( ,Y 



s=0 



Z N 



Ay4+\g ph \*s/(ZN) 
I 



zz(i-z u ) 



N-s 



(59) 



For N large enough the binomial distribution can be ap- 
proximated by a Gaussian 



those known from the standard analysis of collapses and 
revivals of coherent-state evolutions 23] . 



N 



Z10--ZJ) 



N- 



g 2NZ(1-Z) 



y/2iiNZ(l - Z) V2nNZ 



whose shape is controlled mainly by the product NZ. 
The smaller Z the less important its exact value (the 
second approximate equality holds for small Z). So NZ 
is the parameter we may have a chance of seeing in ex- 
periments. Let us stress that both N and Z are relativis- 
tically invariant. 

The atomic inversion w(t) is here a sum of N oscilla- 
tions, each at a different Rabi frequency. It is clear that 
for finite N the evolution of w(t) will reveal collapses 
and revivals, and not a simple Rabi oscillation as would 
be expected on the basis of irreducible representations 
(cf. the analysis of the experiment in the next subsec- 
tion, and in particular Fig. 6). The solution is similar to 



C. Experiment: Thermal mixture as the initial 
condition 



Replacing the Poisson statistics in l|51|l by thermal 
probability we obtain a formula applicable to situations 
where the coherent light is replaced by a thermal mix- 
ture. We thus assume (14811 with 



N 



Tr Pfiold n(n, S ) = P(n)( s ) Z^l - Z^-% (60) 



P(n) 



(61) 



(l + n)(«+!)' 

and find the formula directly applicable to experiments: 



J 



w(t) 



P+ -P- 



N 

E 



N 



z s u (i - z u 



oc 



Taking the limit TV — > ex: we obtain 



P- U \\ 12 

-p+)\g\ 



n=0 



n=0 



(n+ l )s sin 2 (ty/% + \g\*<^) 
N A 2 /4+ \g\ 2 (n + l)s/N 



, |2 sin 2 (t^A 2 /4+| 5ph | 2 (n + l) X , 
P+)\9ph\ X^- 



A 2 /4+|<? ph | 2 (n + l) X . 



.(62) 



(63) 



r 



which, up to Xuj is known from irreducible representa- In what follows we compare theoretical results with the 
tions. precise data on optical Rabi oscillations reported by the 
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FIG. 1: Data from the experiment vs. standard theoretical 
predictions. The blue line is the Rabi oscillation in cavity 
with photon lifetime T cav = 220 /is. The red curve fits better 
but corresponds to T cav = 45 /js. 
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FIG. 4: The same parameters as in Fig. 3 but for longer times 
and in an ideal cavity. 



.4 
0.2 



w(t) 



0.2 0.4 0.6 0.8 



■- wvwWjWW^w^ --^ t[ms] 



1 



-0.2 
-0 .4 



FIG. 5: The same parameters as in Fig. 3 and Fig. 4 but 
additionally damped by exp(— t/T ca ,v), with r cav = 500 fis. 



FIG. 2: Thermal probability distribution with n = 0.05 

Kastler-Brossel Laboratory group from Paris @ . 

The relevant plot is Fig. 2 A in @. In Fig. 1 we show 
the data vs. standard theoretical predictions based on 
irreducible representations. We assume exact resonance 
condition A = 0, g p h — 47 kHz, and p + = 0.99. As the 
initial field state we take the thermal mixture P(n) for 
n = 0.05 (Fig. 2 shows that P(n) agrees with Fig. 2a 

w(t) 




t [fJ-S 



FIG. 3: Data from the experiment vs. predictions based on 
the reducible representation with NZ = 28. The blue curve 
corresponds to an ideal cavity. The relaxation is not a decay 
but a beat: Waiting sufficiently long we will see a revival. 
The red line is the Rabi oscillation with additional damping 
corresponding to T cav = 220 /is. 



from jjj). The blue curve represents the Rabi oscillation 
additionally damped by the factor exp(— t/T cav ), T cav = 
220 /is, which could be expected on the basis of cavity 
relaxation parameters. The red curve is the same Rabi 
oscillation but with stronger damping, T cav = 45 /is. As 
we can see the red curve almost agrees with the data; 
the first two minima are somewhat lower than the data 
but the analysis does not include center-of-mass motion, 
finite duration of experiment, and the resulting frequency 
spread and detunings. 

So the red curve yields a reasonable agreement with 
experiment. The problem is that the cavity lifetime was 
220 (is and not 45 fis, and the data should coincide with 
the blue curve. This observation agrees with the remark 
of Brune et at that 'cavity relaxation plays a marginal 
role in the decrease of oscillation'. The authors further 
write: 'Dark counts (...) are one of the main causes of 
oscillation damping(...). Decoherence by collisions with 
backgroundgas may also contribute to the oscillation re- 
laxation' (|9J, p. 1801). The relaxation observed experi- 
mentally thus appears to be stronger than expected, but 
the authors do not discuss the subtlety in much detail. 
In order to fit the data they used 'damped sinusoids', 
but did not explain what kind of a damping factor (or 
factors?) were employed. Also the Fourier analysis was 
performed on time-symmetrized signals, and not just on 
damped ones. This is important since symmetrized sig- 
nals are closer to revivals-collapses typical of beats than 
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w(t) 




[jJS] 



■t[s] 



-0 . 4 



FIG. 6: The same w(t) at two different time scales, for 
NZ — 1000, <7 p h = 47 kHz, ideal cavity, and ideal vacuum. 
For < t < 100 /ts (duration of the measurement reported 
in the plot is indistinguishable from the ideal undamped 
Rabi oscillation predicted by irreducible representations. For 
< t < 0.05 s we observe collapse and revival. The upper plot 
illustrates the idea behind the correspondence principle: For 
any finite time interval one can choose NZ in a way guar- 
anteeing an agreement, within some given error bars, with 
the standard theory. The lower plot shows that experiments 
involving sufficiently long times are in principle capable of 
discriminating between finite and infinite JVs. 

to damped oscillations. In the light of the results we dis- 
cuss in the present paper it is clear that a more detailed 



For the same reasons as those discussed in the preceding 
subsections we obtain the standard formula but with the 
cut-offs Xuj = Z u /Z, x'uj = t/Xui an d renormalized e p h = 



analysis of experimental damping factors is required. 

In Fig. 3 we show the result with the same mixed initial 
condition but now computed by means of the reducible 
representation. The parameters are A = 0, g p h = 47 
kHz, Z = 0.1, and N = 280. The blue curve involves 
no damping (ideal cavity); the red curve is additionally 
damped by exp(— t/T cav ), T cav = 220 (is. The first two 
minima are again lower than the data, but the general 
agreement with experiment is acceptable. As one can see 
the damping is caused here mainly by the beats. Un- 
fortunately we did not possess all the data needed for a 
realistic comparison with the experiment. 

Now, can we distinguish between true damping and 
beats? In principle yes, but we must wait longer and 
have better cavities. Fig. 4 shows the dynamics of w(t) 
in the reducible representation with exact vacuum initial 
condition, ideal cavity, Z = 0.1, N = 280, A = 0, and 
9 P h = 47 kHz for < t < 0.0015 s. Fig. 5 shows the same 
dynamics but now additionally damped by exp(— t/T cav ), 
?cav = 500 /xs. In principle, in a roughly twice better 
cavity we might see the revival. It must be stressed that 
this revival has no counterpart if field is quantized in irre- 
ducible representation, and vanishes if N — > oo (the limit 
taken with Z fixed). Fig. 6 shows the dynamics of w(t) 
in ideal cavity, exact vacuum, and for N = 10000, the 
remaining parameters being kept as before. This case is 
interesting since up to 100 /is, that is the time available 
in the discussed experiment, the plot is practically indis- 
tinguishable from the result based on irreducible repre- 
sentations, whereas for longer times we observe collapse 
and revival. The greater NZ, the later comes the revival. 

The examples show that the origin of relaxation should 
be carefully reexamined. 

D. General coherent-state initial condition 

Let us now replace vacuum by a general \z) but for 
simplicity keep the initial atomic state to be |+). Then 



(64) 



(65) 



eov~Z, z p h = zyZ. The rescaling z i— ► z p h is exactly 
analogous to l|24|l . up to the presence of the cut-off, a 
negligible modification if Xu ~ 1- 



w(t) = (+\Rz(t)\+) 

= 1 >f y |g P h| 2 (» + 1) s sin 2 (VAV4 + |g ph | 2 (n + l)s/(ZN)) |z^| 2 " _ ]z ^ f N \ 
2 Z N A 2 /4+|5 P h| 2 (n + lMZA0 n\ e \s) ^ 

The limiting form, for TV — > oo, is again familiar 



n^oo 2 ^ A 2 /4+ \g ph \ 2 (n + l)x, 

I 



w(t) 



1 
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i f% A 
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2 
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y 20 ^° w w ' 



FIG. 7: w(t) computed on the basis of irreducible represen- 
tations for a coherent state vs. experiment (Fig. 2B in 9]), 
with 2 = s/QA, T = K, g p h =47 kHz. Cavity param- 
eters: T cav = 50 [is (red), T cav = 220 /is (blue), and ideal 
cavity (green). The two level system is in a mixed state with 
p+ = 0.97. The plot roughly coinciding with the data corre- 
sponds to a cavity that was more than four times worse than 
the one actually used in the experiment. 



w(t) 




t [/US] 



FIG. 8: w(t) (red) computed on the basis of the reducible 
representation, with N = 420, Z = 0.1, and T cav = 220 /is, 
for the same coherent state as in Fig. 7. The blue line is 
the plot based on the standard formalism with T cav = 50 /is. 
The reducible formalism produces the desired result with a 
realistic value of the damping parameter. 



We again compare theory with the experiment of the 
Paris group. Fig. 7 shows the predictions based on the 
standard (irreducible) formalism and corresponding to 
Fig. 2B from We have the same problem as before: 
In order to fit the data we have to assume the cavity 
lifetime T cav = 50 fis, which is much worse than in the 
experiment. On the contrary, the reducible formalism 
leads to a reasonably looking curve even for T cav = 220 
/is if we choose NZ = 42 (Fig. 8). In these data we did 
not have access to the error bars. Fig. 9 and Fig. 10 show 
the plots corresponding to Fig. 2C in [jj. Finally, Fig. 11 
illustrates the correspondence principle: Given any data 
collected in a finite time interval the reducible representa- 
tion is capable of reconstructing the experimental points 
with arbitrary precision. 
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w (t) 




FIG. 9: w(t) computed on the basis of irreducible represen- 
tations for a coherent state vs. experiment (Fig. 2C in 9]), 
with z = V0.85, T = K, g p h = 47 kHz. Cavity parameters: 
(a) T cav = 50 ^s (red), (b) T cav = 220 fxs (blue), and (c) ideal 
cavity (green). The two level system is in a mixed state with 
p+ = 0.99. 
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FIG. 10: w(t) (red) computed on the basis of the reducible 
representation, with N = 420, Z = 0.1, and T cav = 220 /is, 
for the same coherent state as in Fig. 9. The blue line is the 
plot based on the standard formalism with T cav = 50 /j,s. 



VII. CONCLUSIONS 

The main conclusions are the following. First of all, the 
evolution of atomic inversion during a finite time interval 
can always be reconstructed by means of some N < oo 
representation. So, finite N representation can be re- 
garded as a generalization of standard quantum optics. 
This is a consequence of the fact that there exists the 
correspondence principle N — > oo, analogous to h — > or 
c — ► oo. To our surprise at least some part of the avail- 
able data seems to be more consistently explained with 
finite N than with the standard formalism. In order to 
find the proof that N < oo is physical, we need to find a 
revival of a decaying vacuum Rabi oscillation. 
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FIG. 11: Correspondence principle in action. The green curve 
represents prediction of the standard theory for w(t) with the 
coherent state initial condition, z = V0.85, p+ = 1, in ideal 
cavity and zero temperature. The remaining curves are the 
predictions based on reducible representations with Z = 0.1, 
z ph = y/OM, for: TV = 200 (black), TV = 2000 (blue), and 
TV = 10000 (red). 



Appendix A: 'Reducible field quantization' in 
questions and answers 

Here we have collected the points we think are crucial 
for a correct understanding of the formalism based on 
TV < oo representations. 



A. What if Z u = 1? 

In principle this case is not excluded. Then the vacuum 
is monochromatic, i.e. 



\0) = 
10 = 



|p,0), 
|p,0}<8> 



|p,0>, 



(66) 
(67) 



Z = 1, and Xuj = 1 (exactly). In such a vacuum a photon 
with momentum different from p cannot occur. In partic- 
ular no resonance fluorescence is then possible. The law 
of large numbers is then trivial since we are dealing with 
the Bernoulli process with probability of success equal to 
1. The parameter TV then cancels in all the formulas for 
w(t). However, this state is not very physical and there 
is no reason to believe such a vacuum can be encountered 
in experiment. 



The cut-off is a consequence of square-integrability of the 
wave function. Let us add that the wave function 0(k) 
has a status analogous to the wave function of the Uni- 
verse. 



C. How small is Z? 

We do not know. However, if differences between 
x(fe) and 1 are negligible in the quantum optics regime 
then many different momenta are equally probable and 
Z (which is the maximum of the probability) must be 
a very small but nonzero number. Then, on the other 
hand, TV must be very large, and practically the param- 
eter that controls finite TV representations is the product 
NZ which, as we have seen, may be of the order of hun- 
dreds or thousands. In a sense, the smaller Z the better 
since then instead of two parameters, TV and Z, we have 
a single one: NZ. 



D. What is the link between N < oo 
representations and oscillator wave packets? 

Take a nonrelativistic oscillator with the Hamiltonian 



H = P^/2m + mCl A Q 2 /2 



(68) 



and assume that CI is not a parameter but an operator 
which nevertheless commutes with P and Q. This mod- 
ification seems trivial, but is not if CI is nontrivial. If 
CI = loI u is its spectral representation with spectral 
projectors I u , then [I U ,P] = [I U ,Q] = [Iu,H] = 0. De- 
fine Hu, = L,H. Then 



H 

J 



x huj 



= I. 



(69) 



(70) 

(71) 
(72) 



B. What is the meaning of \^ w 1? 

What it means is that the probability Z(k), treated 
as a function k t— > Z(k), is very flat in some part of its 
domain (termed in [y, 0,13 the quantum optics regime). 
Our resonant frequency is assumed to belong to the flat 
region of the distribution, i.e. x(p) is close (or even 
equal) to 1. The function x(fc) = Z(k)/Z has then all 
the features typical of the cut-off functions employed in 
quantum optics. However, Z[k) cannot be a constant 
function since then the wave function normalization con- 
dition J2k z ( k ) = Efc |0(fc)| 2 = 1 would not be fulfilled. 



This is precisely an example of N = 1 representation. 
Taking TV such (noninteracting!) oscillators we arrive at 
a TV < oo representation. The oscillator can exist in 
superposition of different eigenstates of CI. A gas con- 
sisting of TV such oscillators, all in identical ground-state 
wavepackets, is a Bose-Einstein condensate. But this is 
simultaneously our vacuum state. 



E. What about Lorentz invariance? 

No problem. A relativistic formalism is developed in 
I (although is yet preliminary — one should 
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work with Q). We plan to redo the Jaynes-Cummings 
calculations in the representation of [8(. 



F. Are there modifications of the blackbody 
radiation? 

Not really. The analysis of the problem given in Q 
was premature since the role of the large- N limit and 
the law of large numbers was not yet understood in this 
context. A single oscillator, i.e. the N = 1 represen- 
tation, has spectrum identical to the standard one (i.e. 
E n = fojj{n + 1/2); the difference is that here u) is also 
an eigenvalue and not a parameter) . The held is a gas of 
N such oscillators. The only subtlety is that the gas is 
finite and, as shown in [j| , the statistics may be based on 
Renyi a-entropies with a — 1 — 1/N. The limit N — > oo 
is then equivalent to a — > 1, but this is the Shannon limit 
of a-entropies. 

G. What about the divergences? 

This is not yet completely clear, but all the examples 
discussed so far show that the cut-offs occur in the correct 
places. The work on full QED, loop integrals included, is 
in progress. There is one element that is not completely 
controlled yet. Namely, if we work in full space and not in 
a cavity (such as the classical current example discussed 
in p|) then it is natural to start with N = 1 representa- 
tions that involve spectral projectors 1^ corresponding to 
plane waves (i.e. these are not really projectors because 
of Dirac delta normalizations). Such a procedure is useful 
in some cases, but the final formalism should always use 
actual projectors in order to avoid introducing artificial 
infinities coming from the terms 1^. The 'modes' should 
be associated with basis vectors in a Hilbert space, and 
not with plane waves. A trivial way out is to work al- 
ways with finite volumes, but the formalism then lacks 
elegance. Some polishing of the formalism is here yet 
needed. 



H. What about vacuum energy? 

For a finite N the vacuum part of the Hamiltonian is 
a well defined operator from the center of the algebra 
(i.e. commutes with everything). We can remove it by a 
unitary transformation which is well defined (this is what 
we implicitly do in the present paper, and sometimes refer 
to the procedure as a vacuum picture |||). The average 
energy of a single oscillator is finite by assumption (this 
is a condition on the domain of the N = 1 Hamiltonian 
and means that ^ w ujZ^ < oo). For an arbitrary N the 
vacuum energy is N times the average energy of the N — 
1 case and, of course, diverges with N — > oo. However, 
in this sense the mass of a glass of water diverges if one 
treats the thermodynamic limit too literally. By the way, 



the vacuum energy of Dirac electrons, as discussed in |7(, 
is negative. The discussion of vacuum energy in QED 
must involve both fermions and bosons, and then we have 
a difference of two finite expressions which may be well 
defined even in the limit of large N. 



I. Isn't what we do a cut-off regularization in 
disguise? 

No, because there are no cut-offs in operators, and in 
the Hamiltonians in particular. It is true that in effect 
the end result is similar, especially in the limit N — ► oo. 
But since there is no cut-off in the Hamiltonian, there 
cannot be any cut-off dependence in its spectrum! This 
point is very important since, in principle, it can lead to 
yet another direct test of the N < oo representation. In 
a forthcoming paper we shall discuss the spectral prop- 
erties of the Jaynes-Cummings model. The experimen- 
tal aspects have been worked out by the Caltech group 

mm. 

Appendix B: Technicalities 

J. Spectral projectors of a(s)^a(s) vs. a(s) and a(s)t 

For s — we find a(s) = 0. So consider s > 0, to > 0, 
and 

N t 

n(m, s)a(s)Il(n, s) = H(m 7 s)a(sya(s)a(s)Il(n,s) 

ms 

= n(m, s)a(s)H(n, s). 

to 

Hence n -^ m Il{m, s)a(s)JJ(n, s) = and for < to ^ 
71—1 one finds 

n(m, s)a(s)TL(n, s) = 
and, by Hermitian conjugation, 

II(n, s)a(s) + n(m, s) = 0. 
Analogously, for to = 0, n > 0, 

11(0, s)a(s)U(n, s) = -U(Q,s)a(s)Ii(n,s) 
n 

and 

n — 1 

n(0, s)a(s)U(n, s) = 0. 

n 

Now assume a(s)n(0, s) ^ 0. Then there exists a vector 
\tf>) such that a(s)H(0,s)\ip) =: |*) ^ 0. However, 

(*|*) = (ip\Il(0,s)a(s) f a(s)Tl(0,s)\ip) = 0. 

Contradiction. Therefore a(s)H(0, s) = 0. 
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K. Spectral projectors of X(s) vs. Q(s) 
For n > 

[R+a(s),U{n,s)] = R+(a{s)IL(n, s) - II(n - 1, s)a(s)). 

Employing the formulas from the previous subsection we 
find 

a(s)H(n,s) — II(n — 1, s)a(s)H(n, s), 
II(n— 1, s)a(s) — II(n — 1, s)a(s)H(n, s), 

and [R+a(s), II(n, s)] — 0. By Hermitian conjugation 
[.R_a(s)t, ft(n, s)] =0. For n = 

[i? + a(s),fi(0,s)] = i?+a(s)n(0, s) = 0. 

By conjugation [.R_a(s) T , 11(0, s)] = and, for all n > 0, 
[ft(s),I3(n,s)] = 0. 

L. Average (z| II (n, s)| 2) 

By construction a(s)\Q) — 0. Introducing the rescaled 
operators a(s) = a(s)/ y/s/N, satisfying [a(s) ) d(sy] = 



P(s/N), we check that 



a( s ya(s)(a( s y) n \0) = (£) ^(,s)t 5 ( s )(a(,s)t)"|0) 
= (^) 1+f n(a( S )t)"|0) = ^n(a( S )t)"|0). 



Therefore 



n(n', S )(a( S )t)"|0> = 5 nin ,(a(s?) n \0), 

U(n,s)\z) = (s/Ar) (a(*)TlO>, 



and 



(z\U(n,s)\z) = ( s /N) nl -^—e-^^ s/N) (0\P(s/N)\0) 



2 / JV 



N-s 



J 



M. Explicit form of R 3 (t) 



We have to simplify 

R 3 (t) 



(73) 



cos (V AV4 + M^ ) - ; s '" U A" 1 ' U "'^ 



One checks that [A, i? 3 ] = [X,ft] = [A,i?+a w ] 
[X, i?_a* ] = [X, a* a w ] = 0, [ft, i? 3 ] = gR-a* u - gR + a u 



Employing these formulas we arrive at (|15fl . 



fti? 3 ft - (A 2 /4 - \g\ 2 X)R 3 + -(gR+a^ + gR_a* u 
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